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Abstract 

This  paper  develops  local  signal  detection  strategies  for  spectral  resolution  of  frequencies  of 
nearby  tones.  The  problem  of  interest  is  to  decide  whether  a  received  noise-corrupted  and  discrete 
signal  is  a  single-frequency  sinusoid  or  a  double-frequency  sinusoid.  This  paper  presents  an  extension 
to  M.  Shahram  and  P.  Milanfar  (On  the  resolvability  of  sinusoids  with  nearby  frequencies  in  the  pres¬ 
ence  of  noise,  IEEE  Trans.  Signal  Process.,  to  appear,  available  at  http://www.soe.ucsc.edu/~milanfar) 
the  case  where  the  noise  variance  is  unknown.  A  general  signal  model  is  considered  where  the  fre¬ 
quencies,  amplitudes,  phases  and  also  the  level  of  the  noise  variance  is  unknown  to  the  detector.  We 
derive  a  fundamental  trade-off  between  SNR  and  the  minimum  detectable  difference  between  the 
frequencies  of  two  tones,  for  any  desired  decision  error  rate.  We  also  demonstrate  that  the  algorithm, 
when  implemented  in  a  practical  scenario,  yields  significantly  better  performance  compared  to  the 
standard  subspace-based  methods  like  MUSIC.  It  is  also  observed  that  the  performance  for  the  case 
where  the  noise  variance  is  unknown,  is  very  close  to  that  when  the  noise  variance  is  known  to  the 
detector. 
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1.  Introduction  and  problem  setup 

Resolving  sinusoidal  signals  with  nearby  frequencies  has  been  of  a  significant  interest 
in  array  processing  and  in  particular  direction  finding,  where  two  incoherent  plane  waves 
are  incident  upon  a  linear  equi-spaced  array  of  sensors  [2].  The  previous  works  done  in 
this  area  fit  in  two  categories;  some  have  researched  and  developed  novel  methods  and  al¬ 
gorithms  (see  Refs.  [3,4]  for  a  list  of  the  related  literature  and  Refs.  [5-7]  for  more  recent 
works),  whereas  many  others  have  focused  on  the  performance  analysis  of  the  developed 
algorithms  [2,8-16].  Some  common  approaches  in  the  latter  group  have  been  to  derive  a 
sensitivity  measure  for  the  algorithms  to  the  noise  level  [14],  or  to  determine  the  “thresh¬ 
old”  SNR  at  which  a  minimum  resolvability  (in  statistical  terms)  can  be  obtained  [2,15]. 
In  any  event,  the  trade-off  between  resolvability  and  SNR  has  been  consistently  used  as  a 
performance  figure  for  spectral  estimators. 

We  consider  the  signal  of  interest  to  be 

^(x;  (5i,  (52)  =  fli  sin(27r(/c  -  Si)x  +  (pi)  +  a2sin(27r(/c  -f  S2)x  +  (P2)  (1) 

with  X  e  [—B/2,  B/2],  where  the  frequencies  of  sinusoids  (fc  —  (5i  and  fc  +  S2)  are  around 
a  “center”  frequency  /c.  This  center  frequency  can  be  assumed  to  be  known  or  estimated 
beforehand  by  applying  one  of  the  existing  spectral  estimation  methods.  Assuming  that  the 
measured  signal  is  corrupted  by  noise  and  sampled  at  rate  of  /s,  we  can  write  it  as 

f(k-,Si,S2)^s(k-,Si,82)  +  w(k)  (2) 

=  ai  sin|^27r(/c  -  (50^  +  4>i^ 

+  a2sin^7t(fc  +  S2)Y+(p2^  +  w(k),  (3) 

where  the  integer  index  k  is  in  the  range  k  e  {—{N  —  l)/2, . . .  ,{N  —  l)/2}  and  A  =  B/s 
is  the  total  number  of  samples.  The  term  w(k)  is  a  zero-mean  Gaussian  white  noise  process 
with  unknown  variance  cr^. 

The  spectral  representation  (i.e.,  the  discrete-time  Fourier  transform  (DTFT))  of  the 
signal  in  (3)  consists  of  two  overlapping  patterns  centered  at  /c  —  Si  and  /c  -h  S2.  According 
to  the  so-called  Rayleigh  criterion  [10],  these  two  peaks  in  the  frequency  domain  are  barely 
resolvable  if 

Si+S2=1/B.  (4) 

Hence  for  any  signal  with  a  frequency  separation  (Si  -|-  S2  <  l/B),  the  main-lobe  of  the 
Fourier  transform  of  the  (sum  of)  two  sinusoids  is  located  in  the  same  DTFT  bin.  In  other 
words,  the  two  frequency  components  are,  in  the  Rayleigh  sense,  unresolvable .  This  is  the 
scenario  of  interest  in  this  paper.  In  the  remainder  of  this  paper,  we  use  the  phrase  “signals 
with  short  observation  interval”  to  identify  signals  in  which  the  values  of  B,  Si,  and  S2 
satisfy  the  inequality  Si  4-  S2  <  l/B. 

As  we  are  to  decide  whether  the  received  signal  is  a  single-sinusoid  (one  peak  in  the 
spectral  domain)  or  double-sinusoid  (two  peaks  in  the  spectral  domain),  the  problem  of 
resolution  can  be  posed  as  the  following  hypothesis  test 
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Hq:  (5i  =  0  and  82  —  0,  ,,, 

Hi:  3i  >0or52>0, 

where  Ho  and  Hi  denote  the  null  hypothesis  (one  peak  is  present)  and  alternative  hy¬ 
pothesis  (two  peaks  are  present),  respectively.  Since  we  consider  the  case  where  ^i  and  82 
are  unknown  to  the  detector,  (5)  represents  a  composite  (but  one-sided)  hypothesis  testing 
problem  [17].  We  treated  this  problem  in  the  case  where  the  noise  variance  was  known 
in  [1],  and  the  contribution  of  this  paper  consists  of  extending  the  earlier  analysis  to  the 
unknown  noise  variance  case. 

Almost  all  the  previously  developed  techniques  in  this  area  employ  the  structure  of  the 
second  order  statistics  of  the  signal.  The  key  assumption  is  that  of  independent  uniformly 
distributed  phases  of  each  sinusoid.  Assuming  cpi  and  (/>2  to  be  independent  and  uniformly 
distributed  random  variables  in  the  range  of  [0, 27r],  the  resulting  hypothesis  test  from  (5) 
using  the  signal  covariances  is  given  by 

Ho:  f~AA(0,Ro  +  a2l), 

Hi:  f~AA(0,Ri -fa^I), 

where  Ro  and  Ri  are  the  autocorrelation  matrices  of  the  signal  s(k)  in  (2), 

Ro  =  Re[r(fc)P(fc)],  (7) 

2  2 

Ri  =  ^  Re[r(/e  -  80Pif,  -  3i)]  +  ^  Re[r(/c  +  82)r'^(fc  +  32)],  (8) 

where  Re[-]  denotes  the  real  part  and  r(  )  is  the  vector  form’  of 

r(k;  fc)  =  exp|^ j27r/c  .  (9) 

For  the  most  idealistic  case  where  all  the  parameters  in  (7)  and  (8)  are  known  to  the  detec¬ 
tor,  an  NP  detector  for  (6)  decides  Hi  if 

rc(f)  =  f[iRi  +  a^I)-’  -  (Ro  +  a2i)-i]f  >  j,,  (10) 

where  subscript  “c”  denotes  the  “completely  known”  case.  In  practice,  however,  since  only 
the  time  series  are  available,  an  estimate  of  the  autocorrelation  matrix  derived  from  the  sig¬ 
nal  samples  is  used.  Furthermore,  since  finding  the  maximum  likelihood  (ML)  estimation 
of  the  autocorrelation  matrix  is  a  highly  complicated  task,  other  suboptimal  alternatives 
are  used.  The  so-called  subspace  methods  (e.g.,  MUSIC)  for  spectral  estimation  are  based 
on  the  eigen-decomposition  of  the  autocorrelation  matrix  into  orthogonal  signal  and  noise 
subspaces  [2^,18-20]. 

Our  proposed  analysis  is  based  on  the  model  for  the  measured  signal  in  (2)  instead  of 
relying  on  the  second  order  statistic  (i.e.,  covariance  structure)  of  the  signal.  It  is  useful 
to  mention  that  our  methodology  assumes  the  phases  of  sinusoid  to  be  deterministic  un¬ 
known  variables,  whereas  the  subspace  methods  treat  the  phases  of  sinusoids  as  uniformly 
distributed  random  variables. 


*  Superscript  “T”  in  (7)  and  (8)  denotes  conjugate  transpose. 
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Our  approach  is  to  quantify  a  measure  of  resolution  in  statistical  terms  by  addressing  the 
following  question:  “What  is  the  minimum  separation  between  frequencies  of  two  nearby 
tones  (maximum  attainable  resolution)  that  is  detectable  at  a  particular  signal-to-noise  ratio 
(SNR),  and  for  pre-specified  probabilities  of  detection  and  false  alarm  (Pd  and  Pf)?” 

Addressing  the  above  question  will  provide  a  fundamental  performance  bound  which 
can  be  used  to  understand  the  effect  of  SNR  and  also  other  parameters  on  the  achievable 
resolution  in  spectral  analysis.  Furthermore,  the  final  computed  performance  limit  is  simply 
the  result  of  employing  a  local  detector  which  can  be  indeed  used  and  implemented  in 
practice.  For  this  purpose,  we  put  forward  some  comparisons  between  a  practical  setup 
of  the  proposed  algorithm  and  the  MUSIC  algorithm.  We  demonstrate  that  the  proposed 
detectors  yield  noticeably  improved  performance  in  resolving  the  spectra  of  nearby  tones. 

In  our  earlier  work  [1],  we  studied  the  problem  in  the  case  where  the  noise  variance 
is  known  to  the  detector.  In  this  paper,  the  case  of  unknown  noise  variance  is  considered, 
which  is  perhaps  a  more  practical  scenario.  The  main  result  of  the  forgoing  analysis,  as  we 
shall  see,  is  that  there  is  little  loss  in  performance  when  the  noise  variance  is  unknown. 

In  Section  2  we  introduce  our  approach  for  the  most  general  case,  where  the  ampli¬ 
tudes  and  phases  are  unequal  and  unknown  to  the  detector.  In  this  section  we  also  develop 
the  corresponding  detection  strategies  and  characterize  their  performance.  In  Section  3 
we  present  some  results  and  comparisons  of  the  proposed  method  with  existing  subspace 
methods.  Finally,  in  Section  4,  we  summarize  the  results  and  present  some  concluding 
remarks. 


2.  Detection  theoretic  approach 

We  consider  the  most  general  case  of  the  signal  model  of  (3),  with  unknown  amplitudes, 
phases,  and  unknown  frequency  parameters  ((5i  and  Sj)  and  also  unknown  level  of  noise 
variance  (cr^).  The  case  of  known  and  equal  amplitudes  and  phases  has  been  studied  in 
[1]  when  is  known,  the  result  of  which  has  been  shown  to  lead  to  a  uniformly  most 
powerful  test. 

When  two  spectral  peaks  corresponding  to  the  frequencies  of  sinusoids  (fc  —  and 
fc  +  Sj)  are  located  in  the  same  bin  (below  the  Rayleigh  limit),  the  range  of  interest  for  the 
values  of  i5i  and  §2  is  small  (Si,  S2  <  1  /2B).  Hence,  it  is  quite  appropriate  for  the  purposes 
of  our  analysis  to  consider  approximating  the  model  of  the  signal  around  (Si,  S2)  =  (0,  0). 
The  second-order  Taylor  expansion  of  (2)  about  (Si,  S2)  =  (0,  0),  with  all  other  variables 
fixed,  is 

2 

s(t,Si,S2)^J2‘^iPi  (k)  +  Pi  qt  (k),  (11) 

i=0 

where 

Piik)^(^Y^  sin^Ttfcy^,  (12) 

„«:)  =  (4)  cos(2,a|-), 


(13) 
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and 

(14) 

(15) 

We  elect  to  keep  terms  up  to  order  2  of  the  above  Taylor  expansion.  This  gives  a  more 
accurate  representation  of  the  signal  since  in  some  cases  the  first-order  terms  (pi{k)  and 
qi(k))  would  be  very  small  or  even  would  simply  vanish.  Rewriting  (11)  in  vector  form 
will  result  in 
2 

-f  A  qi,  (16) 

i=0 

where 

P/  =  [pi{-(N-l)/2),...,pi{(N-l)/2)f,  (17) 

q/  =  [qi{-{N-l)/2),...,qiaN-l)/2)Y.  (18) 

Now,  the  hypotheses  in  (5)  appear  in  the  following  form; 

'  Hq:  z  =  aopo  -f  /6oqo  +  w, 

2 

Tfi;  z  =  -f  ;S,q; -f  w, 

!=0 

where  z  denotes  the  approximate  measured  signal  model.  Equation  (19)  leads  to  a  linear 
model  for  testing  the  parameter  set  6  defined  as  follows: 


(2jty 


i\ 

(2ny 


ai&\  cosj  (pi  +i  — 


-|-fl252‘^os(  (pi  +  i  — 


IT 


ll 


n 


fliSj  sinj  (pi+i  — 


-I- £i2^2si>^(  </’!  +  1^ 


n 


z  =  -1-  w. 

(20) 

H  =  [po  1  qo  1  Pi  1  qi  1  P2  1  q2], 

(21) 

9  =  [ao  Po  ai  Pi  “2  Pif , 

(22) 

where  H  and  ^  are  an  N  x  6  matrix,  and  a  6  x  1  vector,  respectively.  The  corresponding 
hypotheses  are 


Ho:  A0  =  O,  cr2>o, 

Hi:  Ae^O,  a^>0, 

where 

-0  0  1  0  0  0- 

0  0  0  1  0  0 

^“000010 

_0  0  0  0  0  1_ 


(23) 


(24) 


The  hypothesis  test  in  (23)  is  a  problem  of  detecting  a  deterministic  signal  with  unknown 
parameters  (0  and  cr^).  The  generalized  likelihood  ratio  test  (GLRT)  is  a  well-known  ap¬ 
proach  to  solving  this  type  of  “composite”  hypothesis  testing  problem  [21].  The  GLRT 
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uses  the  maximum  likelihood  (ML)  estimates  of  the  unknown  parameters  to  form  the  stan¬ 
dard  Neyman-Pearson  (NP)  likelihood  ratio  detector.  The  GLRT  for  (23)  [21]  gives  the 
following  test  statistic: 

_  ?V[A(HTH)-1aT]-1a? 

zT[I-H(HTH)-1HT]z  ^  ^ 

where  I  is  the  identity  matrix  and^ 

?=(H^H)“*H^z  (26) 

is  the  unconstrained  maximum  likelihood  estimation  of  0.  For  any  given  data  set  z,  we 
decide  Hi  if  the  statistic  exceeds  a  specified  threshold, 

T(z)  >  y.  (27) 

The  choice  of  y  is  motivated  by  the  level  of  tolerable  false  alarm  (or  false-positive)  in  a 
given  problem,  but  is  typically  kept  very  low.  From  (25),  the  performance  of  this  detector 
is  characterized  by  [21] 


(28) 

(29) 

A  =  ^6»W[A(H’^H)-*A'r]“‘A0,  (30) 

where  Q  is  the  right  tail  probability  for  a  central  F  distribution  with  4  numerator 
degrees  of  freedom  and  N  —  6  denominator  degrees  of  freedom,  and  Qpr  is  the  right 
tail  probability  for  a  non-central  F  distribution  with  4  numerator  degrees  of  freedom  and 
N  —  6  denominator  degrees  of  freedom,  and  non-centrality  parameter  X.  For  a  desired  Pd 
and  Pf,  we  can  compute  the  required  value  for  the  non-centrality  parameter  from  (28) 
and  (29).  We  call  this  value  of  the  non-centrality  parameter  A(Pf,  Pd)  as  a  function  of 
desired  probability  of  detection  and  false  alarm  rate.  This  notation  is  key  in  illuminating  a 
very  useful  relationship  between  the  SNR  and  the  smallest  separation  (i.e.,  (^i,  S2))  which 
can  be  detected  with  high  probability,  and  low  false  alarm  rate.  From  (30)  we  can  write 

0-2= - i - 0TA'r[A(HTH)“iA'^l“‘A6».  (31) 

A(Pf,Pd)  L  ^  ^  J 

Also,  by  defining  the  output  SNR  as 


SNR  =  (32) 

(7^ 

and  replacing  the  value  of  cr^  with  the  right  hand  side  of  (31)  the  relation  between  the 
parameter  set  0  and  the  required  SNR  can  be  made  explicit. 


SNR  =  A(Pf,  Pd) 


6>TaT[A(HTH)-iaT]-iA6>  ' 


(33) 


^  Note  that  (HTh)-IhT  is  the  pseudo  inverse  of  H. 
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It  is  instructive  to  simplify  (33)  by  approximating  the  elements  of  the  matrix  H^H.  These 
approximations  yield  (see  Appendix  A  for  details) 


0 


0 


0 


N 

~1 


0 


-A 

0 


0 


0 


-N 

0 


0 

-N^ 


-N 

0  0 

0  0 

-N^ 

24f^ 

-N^ 

160/4 
0  0 


0 

-N^ 


0 


0 

160/4  - 


where  fj,  —  cos(^^A)/sin(^^).  To  gain  further  insight,  we  consider  a  special  case  by 
assuming  ^  a2&2,  which  results  from  a  proper  choice  of  the  center  frequency  fc  (see 
Ref.  [1]),  and  simultaneously  considering  the  case  where  the  value  of  /i  is  close  to  that 
of  (p2-  After  some  algebra,  replacing  N  /  f^hy  B  and  neglecting  non-dominant  terms,  for 
small  (5i  and  &2  ^i,  32  1  /B)  (33)  will  reduce  to: 


45  A(Pf,  Pd)  (ai+a2)^ 

7r4  p4  {ax&\  +  a2&'l)'^ 


(34) 


The  required  SNR  is  minimized  when  3i  =  ^2  =  3  (i.e.,  ai—a2—  1);  and  in  this  case  we 
have 


45  XjPf,  Pd) 
4  (P3)4 


(35) 


An  important  question  is  to  consider  how  different  this  obtained  performance  is  from  that 
of  the  ideal  detector,  to  which  all  the  parameters  (amplitudes,  phases,  frequencies,  and 
noise  variance)  are  known.  We  first  note  that  in  this  case  the  hypothesis  test  in  (19)  will 
be  a  standard  Gauss-Gauss  detection  problem.  Also,  we  can  further  simplify  the  problem 
by  seeing  that  the  term  aoPo  +  /oQo  is  a  common  known  term  under  both  hypotheses  and 
can  be  removed.  As  a  result,  the  following  relationship  can  be  verified  for  the  completely 
known  case: 


SNRid  =  >?(Pf,  Pd) 


6»'^HTH6> 

0TaTAHTHATA6»  ’ 


(36) 


where  the  subscript  “id”  denotes  the  ideal  case  and  rj{Pi,  Pd)  is  the  required  deflection 
coefficient  [21]  computed  as 


'7  =  (e“‘(^f)-e“‘(^d))", 


(37) 
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where  Q  VO  is  the  inverse  of  the  right-tail  probability  function  for  a  standard  Gaussian 
random  variable  (zero  mean  and  unit  variance).  Comparing  the  expression  in  (33)  to  that 
of  (37), 

SNR  _  A.(Pf,  Pd)  0’^ATahThaTa^ 

SNRid  “  r](P{,  Pd)  0'rAT[A(HTH)-iAT]-iA0 

we  notice  that  ri{P{,  Pd)  <  A(Pf,  Pd)  provided  Pd  >  Pf.  Also,  It  can  be  proved  that 
AH^HAT  -  [A(HTH)-1aT]-i  is  a  positive  definite  matrix.  As  a  result,  as  expected, 
SNR  >  SNRid  always. 


3.  Numerical  results  and  comparisons 

We  first  compare  the  performance  of  the  proposed  detector  for  the  unknown  cr^  case 
with  those  of  the  detector  for  the  known  cr^  case  [1]  and  the  ideal  detector  (36).  Fig¬ 
ure  1  presents  the  performance  figures  of  these  detectors  for  the  case  where  ai  =  a2  and 
Si  =  §2  =  S  (curves  of  2SB  versus  required  SNR).  We  observe  that  knowledge  of  the  noise 
variance  makes  little  difference  to  the  performance  (around  1  dB  in  required  SNR).  It  is 
worth  mentioning  that  the  estimate  of  the  noise  variance  used  in  (25)  is  known  to  be  unbi- 


Pd=0.99  Pfa=10-2 


Fig.  1.  28B  vs  required  SNR  for  known  and  unknown  noise  variance,  a[  =  a2  and  51=^2  =  5. 
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ased  [21].  Comparing  the  ideal  (unrealizable)  detector,  the  GLR  detector  in  (25)  requires 
3-5  dB  more  SNR  to  achieve  the  same  resolvability. 

As  for  the  comparison  to  existing  methods,  we  first  note  that  for  subspace  detectors, 
the  phase  is  typically  assumed  to  be  a  uniformly  distributed  random  variable  in  [0,  27r]. 
However  the  “required  SNR”  computed  in  Section  2  is  in  general  a  function  of  the  phases  of 
the  sinusoids.  Hence  in  order  to  compare  the  obtained  results  to  that  of  subspace  methods, 
we  carry  out  the  following  averaging  for  the  required  SNR  over  the  possible  range  of  <f)\ 
and  (p2 : 


Itt  In 

SNRavg  =  ^  y  J  SNRd^i  d<P2,  (39) 

0  0 

where  subscript  “avg”  denotes  the  averaged  value  and  the  integrand  (SNR)  is  the  right-hand 
side  of  (33). 

Next,  we  simulated  the  behavior  of  the  MUSIC  algorithm  for  resolving  sinusoids  with 
nearby  frequencies.  In  simulation  of  MUSIC,  the  signal  is  declared  to  be  resolvable  if  the 
output  of  MUSIC  produces  two  distinct  peaks  within  an  interval  around  the  true  frequen¬ 
cies  (/c  ±  8).  The  simulations  for  MUSIC  are  carried  out  for  cases  in  which  either  a  single 


Pd=0.99  Pfa=10”^ 


Fig.  2.  28 B  vs  required  output  SNR  for  the  MUSIC  algorithm  and  the  proposed  detector,  unknown  noise  variance, 
ai  =  a2  and  =  ^2  =  5. 
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snapshot,  or  multiple  snapshots,  are  available.  Naturally,  we  consider  the  output  SNR  in 
the  latter  case  as  the  sum  of  SNR’s  of  each  snapshot. 

We  develop  two  different  comparison  procedures.  First,  we  compare  the  performance 
of  MUSIC  against  the  detector  in  (25),  where  we  assume  that  the  center  frequency  /c,  at 
which  we  perform  the  hypothesis  test,  is  known  a  priori.  Since  this  might  be  seen  as  an 
unfair  comparison,  we  have  put  forward  an  alternative  (perhaps  more  practical)  scenario, 
too.  In  this  scenario,  we  hrst  seek  assistance  from  MUSIC  to  estimate  the  center  frequency 
(/c)  and  then  apply  the  proposed  detector  in  (25)  centered  at  the  peak  estimated  by  MUSIC. 

The  results  of  these  experiments  are  shown  in  Fig.  2.  First,  we  observe  that  the  pro¬ 
posed  detector  signihcantly  outperforms  MUSIC  in  both  cases  (using  known  or  estimated 
center  frequency).  More  interestingly,  we  see  that  the  result  of  the  proposed  detector  with 
estimated  center  frequency  (provided  by  MUSIC)  is  very  close  to  the  performance  of  the 
same  detector  with  known  center  frequency.  This  implies  that  the  MUSIC  algorithm  does 
a  very  promising  job  in  locating  the  center  frequency  (i.e.,  the  candidate  location  where  we 
can  perform  a  rehnement  step  using  our  proposed  approach).  Intuitively,  the  reason  for  this 
behavior  is  that  for  the  case  where  a  high  probability  of  resolution  (0.99)  is  considered,  a 
fairly  high  value  of  SNR  should  be  provided.  This  value  of  SNR  will  effectively  guarantee 
a  condition  under  which  the  MUSIC  algorithm  will  produce  the  peak  in  its  spectrum  within 
the  range  of  [fc  —  S,  fc  +  5].  This  observation  is  essentially  in  agreement  with  what  has 
been  observed  in  the  past  about  the  stability  of  MUSIC  for  single-sinusoid  signals. 


4.  Conclusion 

The  problem  of  interest  in  this  paper  has  been  to  establish  a  statistical  analysis  of  attain¬ 
able  spectral  resolution  and  to  propose  its  associated  detection  algorithm.  We  formulate  the 
problem  as  a  hypothesis  test,  the  aim  of  which  is  to  distinguish  whether  the  received  signal 
contains  a  single-tone  or  double-tone.  We  have  considered  the  most  general  case  where  the 
amplitudes,  frequencies  and  phases  of  sinusoids  and  also  the  value  of  noise  variance  is  un¬ 
known  to  the  detector.  This  paper  is  different  and  more  general  than  [1]  in  that  it  assumes 

is  unknown. 

By  utilizing  a  quadratic  approximation,  we  in  fact  carried  out  the  analysis  in  the  con¬ 
text  of  locally  optimal  detectors,  and  developed  corresponding  detection  strategies.  The 
performance  hgure  of  resolution  has  been  quantihed  by  the  following  practical  question: 
“What  is  the  minimum  detectable  frequency  difference  between  two  sinusoids  at  a  given 
signal-to-noise  ratio?” 

We  clearly  observed  that  noise  variance  being  unknown  has  little  effect  on  the  detection 
performance.  This  is  a  useful  observation,  since  in  practice  the  variance  is  often  unknown 
to  the  receiver. 

Also,  the  proposed  locally  optimal  detectors  yield  signihcantly  improved  detection 
of  very  nearby  frequencies,  as  compared  to  the  existing  subspace  methods.  In  terms  of 
implementing  the  suggested  detection  algorithm,  we  merely  need  to  estimate  the  center 
frequency.  Fortunately,  as  we  conhrmed  by  some  experiments,  this  task  can  be  very  effec¬ 
tively  performed  by  using  one  of  the  myriad  of  existing  methods  for  spectral  estimation 
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and  then  running  the  proposed  detector  at  the  estimated  peak.  The  performance  of  such  a 
detector  is  nearly  identical  to  that  of  the  detector  with  a  known  center  frequency. 


Appendix  A.  Computing  the  energy  terms 


In  this  appendix,  we  explain  the  general  process  for  the  approximate  computation  of  the 
energy  terms.  We  will  utilize  the  following  identities  for  the  calculation; 


L  1  - 

=  _ _ 

^  I  —X 

k=0 

^  P  om 

^  ^  dx^ 

^=0  m=l 


1 

1  —  X 

1  9 

2  dx 


(A.l) 

(A.2) 

(A.3) 


Y  sinixk)  cos(xk)  =  -  kP  sin^(xk). 

k  k 

Instead  of  showing  all  the  calculations,  for  the  sake  of  brevity  we  discuss,  as  an  example, 
the  calculation  of  the  term  hQho; 

(7V-l)/2 


hJho  =  4 


E  YY) 

k=-{N-l)l2  \  is  / 

(iV-l)/2 

=  E 

k=-{N-\)l2 

(  4^/c 

=  E  2-exp(7— k 


expl 

/s 


k=-{N-\)l2 


■expl  -J—^k 

/s 

.47r/c 


-l2 


■  exp  -j 


/s 


l_expO'?f^(A+l))  l-exp(-j^f£(A+l)) 

=  2A  -  2 - - 2 - - +  2 


;  2nfc 


1  -exp(y^^) 


1  -exp(-7^) 


=  2A  +  2  -  2 


A  '  -  --FV  y 

1  -  cos(^^)  +  cos(^(A  -  1))  -  cos(  Y^(N  +  1)) 

/s  /s  /s 


1  -cos(^^) 

/s 


sin(?^A) 


=  2A-2 


sin(?^) 

Is 


(A.4) 


Since  sin(x)  ^  1  —  |  —  1 1  for  0  ^  x  ^  tt,  and  <  jr,  by  upper  and  lower  bounding 

the  numerator  and  the  denominator  of  |  C  | ,  respectively,  we  have 


|C|=2 


sin(^A) 

Is 


sin(^) 


(A.5) 
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Thus  for  the  range  of  e  ^  2fc/  fs  ^  1  —  e  (representing  the  range  of  fs  from  just  above  the 
Nyquist  rate  (2/c)  to  1  /e  times  the  Nyquist  rate),  we  will  have  \C\  <  1  /e  and  therefore  for 
€  <5/N 

hj^ho  ^  2N.  (A.6) 

A  similar  approach  can  be  followed  to  compute  other  energy  terms. 
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